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PHYSICAL AND DYNAMICAL STUDIES OF METEORS 


Special Progress Report 


1. SUMMARY 


1. 1 Aims 


This research focused on the interplanetary distribution of radar meteors, freed 
from observational bias. We sought both the general average distribution and the 
distribution of meteor streams with their comet and asteroid associations. 

1. 2 Data 


The Havana, Illinois, "synoptic-year" radar observations and the Havana-Sidell 
simultaneous radar-television observations constituted the data. In addition to being 
the largest and most accurate collection, these data have two unique and crucial advan- 
tages: 

A. There is an adequate sample of very slow meteors, which recombination has 
eliminated from all other surveys. 

B. There is a reliable calibration of the meteor masses. 

The only known selection effect not yet adequately assessed is a possible bias against 
fragmenting meteors. 

L 3 Orbital Distribution 


Low-velocity meteors predominate; the velocity mode in the atmosphere is under 
15 km sec -1 , and the mode in space at 1 a.u. from the sun is under 20 km sec 1 . 
Radiants cluster to the sunward and antisun directions. Moderate eccentricities and 
low inclinations predominate; retrograde orbits are very rare. 
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1.4 Space Density 


We extrapolated the orbital distribution observed at the earth to include unobserv- 
able orbits and used the resultant value to correct the space density of observed orbits. 
The result is meaningful between about 0. 4 and 2. 0 a. u. , and uncertain but suggestive 
farther away. Meteors are concentrated to a relatively thin layer centered on the 
ecliptic plane. In that plane, there is a minimum of density at about 0.7 a. u. , and a 
maximum beyond 2 a. u. 

1.5 Collisions 


If this space-density distribution is in equilibrium, then most meteors must end 
their lives in collisions with each other. Moreover, a simplified theoretical model 
incorporating collisions and the Poynting-Robertson effect agrees with observation. 
Nonetheless, collisions need further study. 

1 . 6 Streams 


With the use of the most powerful stream-search technique in existence, we found 
200 new streams in the synoptic-year data and confirmed the existence of two-thirds of 
the 83 streams detected previously in the 1961—1965 sample of radio meteors. Some 
of the new streams have most uncommon orbits. A new, rich stream with a revolution 
period of more than 30 years has been discovered. Streams of low inclination are often 
detected at both nodes. 

1. 7 Associations with Comets and Asteroids 


A number of known comet— meteor associations are confirmed, and a few new 
possible associations detected. The previously detected associations with the minor 
planet Adonis are strongly reinforced (including a major, broad stream), and several 
possible associations with other asteroids, meteorite Pribram, and a Prairie Network 
fireball are suggested. 
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1. 8 Height- Velocity Diagram 


Plots of both the individual radio meteors and the radio streams fail to exhibit 
the discrete-level structure known to exist for photographic meteors. Some features 
of the height-velocity diagram are, however, detected. 

1. 9 Space Density in Streams 

The mean space density in streams is much below the sporadic density, but the 
central density may significantly exceed the sporadic density. The absolute stream- 
density values based on the statistical model of meteor streams and on the sporadic 
density estimated in this report are in an order-of- magnitude agreement with densities 
estimated from the cometary production rates of solid material. 
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2. INTRODUCTION 


The Radio Meteor Project, initiated at Harvard College Observatory under the 
direction of Professor Fred L. Whipple, and transferred to the Smithsonian 
Astrophysical Observatory (SAO) in 1966, constitutes part of a broad program of 
meteor and meteorite research at the two observatories. Under NASA support via 
contract NSR 09-015-033, the project included research activities in Cambridge, 
Massachusetts, and data collection by an eight-station radar system near Havana, 
Illinois, and by an image -orthicon system at Sidell, Illinois. 

Research under the current contract NAS 1-11204 is directed toward distributions 
of radar meteor orbits in interplanetary space. The specific tasks are the following: 

A. Use of the radar /image-orthicon data to improve assessment of radar obser- 
vational biases. 

B. Determination of meteor orbital distributions from the synoptic-year sample, 
including corrections for observational biases. 

C. Determination of orbital parameters of meteor streams and the density dis- 
tribution of meteors in the streams. 

D. Investigation of associations of meteor streams with comets and with asteroids. 

This report contains results for each task and also unanticipated results on the 
space density of meteors and on collisions between meteors. 
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3. OBSERVATIONS 


3. 1 Radar Meteors 


Multistation radar observations of meteors contain the bulk of the orbital data 
on meteors, nearly equaling photographic observations in accuracy while far out- 
numbering them. (Reliable visual or satellite orbits are still rare. ) The Havana 
system, with more receivers and greater sensitivity than any other, acquired the 
largest and most reliable collection of radar meteor data. We are here directly con- 
cerned with the orbital data but must make extensive use of the data on the meteor's 
interaction with the atmosphere to evaluate selection effects. 

3.2 Synoptic Year 

The primary data used for this report are the radar meteor observations collected 
in Havana, Illinois, during the period known as the synoptic year. The equipment and 
observations are described in the Final Report of contract NSR 09-015-033. 

The Havana facility was extensively refurbished in 1967 and 1968: Measurements 
were made of the antenna radiation patterns, and equipment was installed to calibrate 
transmitted and received signals. The synoptic year comprises nearly uniform cover- 
age of all hours and alternate weeks from December 1968 to December 1969, with the 
Havana system in its best condition. 

A satisfactory observation made with the Havana system yields the meteor's 
radiant, velocity, and deceleration; the location of its trajectory in the atmosphere 
and the distribution of electrons along the trajectory; electron-diffusion rates in the 
atmosphere; atmospheric-wind components; and, possibly, electron recombination 
rates or indications of meteor fragmentation. The meteor's orbit around the sun is 
deduced from the radiant and the velocity, and its mass, from the velocity and the 
electron distribution. The computer facilities at Langley Research Center reduced 
the observations. 
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3.3 Simultaneous Radar-Television Observations 


Beginning in 1969, an image-orthicon system on loan from U. S. Naval Research 
Laboratory was sited at Sidell, Illinois, to observe optically some of the same meteors 
observed by radar from Havana. These observations are described in the Final 
Report of contract NSR 09-015-033, and also by Cook et al . (1972). Because we were 
able to use the much better known masses of optical meteors, these simultaneous 
observations are the best data for calibration of the masses of the radar meteors. 
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4. SELECTION EFFECTS 


4. 1 Radar Biases 


The meteors we observe on radar are not, we know, a random sample of the 
influx to the earth; corrections for the many biases have occupied meteor astronomers 
for years. We have lately made several advances, large and small, in measuring 
these biases. All the corrections used are described in the present section. 

4.2 Observational Geometry 

Elford (1964) determined the probability that the Havana radar would observe a 
meteor as a function of radiant altitude and azimuth. His table of the response func- 
tion of the six stations transmitting on the double trough and receiving on a single 
trough at site 3 and a double trough at other sites was not included in his 1964 report 
but is very similar to the tables that were. We have corrected it for the subsequent 
antenna calibration, as follows. The important differences between the measured 
antenna pattern and the design pattern that Elford used are a 6° decrease in elevation 
of the main beam and an increase in sidelobe sensitivity, in the measured pattern. 
Since sidelobe meteors are routinely rejected from the reductions, they are not con- 
sidered here. We first revised Elford' s radiant- sensitivity contours to delete his 
assumption that the electron line density is proportional to cos Z R , where Z R is the 
radiant zenith distance. Next, we rotated the contours by 6° so as to center them on 
the measured main beam. Then we restored a proportionality between electron line 
density and (cos Z^) * 1 which is more nearly correct for fragmenting meteors. 

The weight given each meteor as a function of radiant azimuth and altitude is the 
reciprocal of the revised radiant sensitivity, except that radiants with sensitivities 
less than 10% of the maximum are omitted altogether. A further weight was then 
applied as a function of radiant declination; this is inversely proportional to the range 
of right ascension within the 10% contour. Since observations were uniformly dis- 
tributed over the day, all radiants observable from Havana are now equitably 
represented. 
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4. 3 Ionizing Probability 


Cook et al . 's (1972) ionizing probability, derived from the simultaneous radar- 
television observations at Havana and Sidell, is adopted here. 

Future analysis of more of the simultaneous observations may revise the ionizing 
probability for fast meteors, which has perforce been extrapolated from lower velocities. 
Such an analysis could not, however, remove the predominance of low-velocity meteors. 

4.4 Diffusion 


The upper bound to meteor heights observed at Havana is set by ambipolar diffusion 
of the ionized column (Southworth, 1973). Thus, some of the faster meteors escape 
observation because they are too high. Figure 4-1 shows the distribution in height for 
various velocity groups among 11, 061 meteors from the synoptic year; it is evident 
that the distributions are truncated at the right for velocities above roughly 40 km sec" . 
Fitting a mean height distribution to the left sides of these truncated distributions, we 
found the fraction of meteors missed; we then derived the following empirical correc- 
tion formula. Meteors with velocities v over 34 km sec * have their weights multi- 
plied by the additional factor 

W D = exp [6.56 (log 10 v- 1. 53) 2 ] . (4-1) 

This correction reaches its maximum = 2. 1 at v = 72. 

4. 5 Initial Radius 


Although the initial radius of the electron column (actually the radius reached in 
a few milliseconds) has appeared significant to several observers of relatively bright 
radar meteors, it does not seem important to the synoptic-year observations 
(Southworth, 1973). We propose that the radii measured with two-frequency or three- 
frequency radars (Baggaley, 1970; Moisya et al . , 1970; Kolomiets, 1971) should be 
understood as resulting from a combination of three processes: (1) diffusion of 
individual ions and electrons into the undisturbed atmosphere, as treated theoretically 
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by Manning (1958) and 5pik (1955); (2) expansion of a hot cloud of neutral gas; and (3) 

lateral spread of fragments of the original meteoroid. Only the first of these has 

usually been considered, but the Havana-Sidell observations appear to show that (3) 

is quite significant for some meteors. Baggaley's conclusion that initial radius is 

proportional to the 0. 5 power of atmospheric mean free path is inexplicable by use of 

(1) alone, but easily understood as (2) or as a combination of (1) and (3) or of (1), (2), 

and (3). For an order-of-magnitude estimate of (2), we compute the radius r of the 

c 

cylinder in the atmosphere that contains the same amount of heat energy as the energy 
deposited in the atmosphere by the meteors. We find 


r 2 »6Xl(T 14 2i. 
c pt 


where m and v are the meteor's mass and velocity, respectively; p and t are atmos- 
pheric density and temperature, respectively; and all units are cgs. At 100-km height 
and 60 km sec the brightest synoptic -year meteors (~10^ elec cm -1 , mass ~ 10 - ^ g 

at 60 km sec 1 ) have r ~ 14 cm, which is negligible at the Havana wavelength (733 cm). 

C 12 
However, at the transition between underdense and overdense trials (line density ~ 10 

elec cm ^), a meteor of the same height and velocity would have r fi - 140 cm and an 


echo attenuation 


A”exp[-8^ 2 0f) 2 ]=O.O6 


which is a severe observational bias. It is probable that process (2) has been neglected 
heretofore because ionizing efficiencies were believed to be higher, and masses corre- 
spondingly lower. 


In view of the foregoing, we do not believe that the initial radius biases the 
synoptic-year data, except perhaps as the effect of fragmentation, which is separately 
discussed. No correction was applied. 

4. 6 Recombination 

Southworth (1973) showed that recombination is an important limitation to surveys 
using less-sensitive equipment. F or the synoptic year, however, recombination can 
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be expected to affect only the lowest bright meteors observed, and then not to attenuate 
them beyond observability. An empirical check comparing the height distributions for 
different ranges of radar magnitude showed no apparent effects of recombination. No 
correction was applied. 

Future studies that treat individual masses will need recombination corrections to 
electron line densities and therefore to masses on at least some meteors. However, 
as explained in Sections 4. 8 and 5. 1, the electron line densities do not enter the weights 
in this report. 

4. 7 Fragmentation 

The Havana-Sidell simultaneous observations (Cook et al . , 1972) appear to show 
that fragmentation is sometimes associated with severe attenuations of radar echo 
strength. This can be easily understood as the effect of lateral spread of the particles 
from the mean trajectory. Simple theory leads to the expectation that the effect is 
most important for slow meteors, and photographic meteor studies lead to the expecta- 
tion that it is important for some showers. However, there is as yet no radar analysis 
evaluating its importance, and no corrections have been applied. 

The possible underrepresentation of fragmenting meteors is the only significant 
bias that we suspect to remain in the results reported here. 

4. 8 Mass Distribution 

3 

Although we have observed meteors covering a range of more than 10 in electron 
line density, the steep dependence of ionizing probability on velocity ensures a strong 
negative correlation between mass and velocity in the synoptic-year meteors. As it 
is not possible to select an unbiased sample from among these, we have assigned to 
each meteor a weight representing the number of meteors of some standard mass. 

This procedure implicitly assumes that the mass distribution is independent of orbital 
parameters over the range of masses observed at any one velocity, and we have not 
noticed any counterexamples in our data. Our data are not sufficient to show whether 
the orbital distribution is independent of mass over the whole range of masses observed. 
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However, we expect that there is some dependence of orbital distribution on mass, and 
our results necessarily refer to the orbital characteristics of masses that are larger 
at low velocity than they are at high. 

On the assumption that orbital elements are independent of mass, a small mass 
is as good a sample of orbital elements as is a large one. Since, as just explained, we 
make this assumption for all meteors of the same velocity, we are also taking advan- 
tage of it to equalize, for all meteors of the same velocity, that part of the statistical 
weight that depends on mass. This amounts to omitting the observed total number of 
electrons from the weights. It has appreciably smoothed the final distributions. 

We have adopted Whipple's (1967) mass distribution, specifically that the number 
of particles with mass ^mis 

n = constant (m)"’*’ ^ . (4-4) 

The corresponding factor in the statistical weights to reduce counts to observations 
at equal masses is 

f=P" 1,36 , (4-5) 

where (3 is the ionizing probability. There is certainly a suggestion in Whipple's 

collection of data that the slope is lower for our smaller masses. Adopting an 

-5 

exponent of -1.0 for masses below 10 g would increase the relative abundance of high- 
velocity orbits in our distributions by a factor ~ 3, but would still leave them numer- 
ically unimportant. Since the ionizing efficiency for high velocities also remains 
uncertain, we have not made this adjustment. 
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5. ORBITAL DISTRIBUTIONS 


5. 1 Synoptic -Year Sample 

This section presents tables and graphs of the velocities, radiants, and orbits of 
meteors observed in the synoptic year. All graphs and tables are in three forms: 

(1) unweighted, (2) reduced to equal masses in the atmosphere, and (3) reduced to 
equal masses in space at the earth's distance from the sun. We use three forms to 
facilitate comparison with the great variety of other published forms. Observations 
from earth satellites should be compared with the second form; and observations not 
restricted to the ecliptic (e. g. , comets), with the third form. 

The tape of reduced observations contains 19, 698 meteors. We omit 505 with 
slant ranges over 400 km because such observations are much more likely than others 
to be confused by wind shears or plasma resonance. Accordingly, the unweighted 
sample is 19, 193. For the sample reduced to equal masses in the atmosphere, we 
omit a further 4973, which were observed outside the main lobe of the antenna pattern 
and which would therefore introduce large statistical noise in the results if they were 
given the very large and uncertain weights nominally belonging to them. Thus the 
atmospheric sample is 14, 220. For the sample reduced to equal masses in space at 
1 a.u. , we omit 1622 more meteors with radiant latitudes south of the ecliptic, since 
we cannot get a complete sample of all southern radiants. This leaves 12, 598 in the 
space sample. 

Because we assume that in our observational dynamic range the orbital parameters 
are relatively independent of electron line density, the weights do not include the 
integrated number of electrons observed, although they do include all other mass 
factors; this procedure has appreciably reduced statistical noise in the results. Since 
the equipment parameters were nearly constant over the synoptic year, there are no 
corrections dependent on date of observation. All other corrections described in 
Section 4 are included. The space sample includes corrections for the collision 
probability with the earth. 
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5. 2 Velocity 


Figure 5-1 shows the distribution of no-atmosphere velocities: dotted lines for the 
unweighted sample, dashed for the atmospheric sample, and solid for the space sample. 
By comparison with distributions observed on less- sensitive equipment (e.g. , the 
earlier Havana data; Hawkins, 1963), there is a considerable enhancement in the fre- 
quency of velocities under 30 km sec*. Moreover, velocities down to 11 km sec * 
are observed, whereas earlier surveys barely reached 15 km sec *. Recombination 
caused the bias against slow meteors on less-sensitive equipment (Southworth, 1973). 

Corrected for ionizing efficiency and reduced to equal masses, the velocity dis- 
tribution in the atmospheric sample shows the greatest concentration toward low 
velocities of any published distribution known to us. The mean is approximately 
14. 5 km sec *. However, the space sample shows that these very low velocities are 
less common outside the earth's gravitational field. 

5. 3 Radiants 


Figures 5-2, 5-3, and 5-4 show the distribution of radiants with respect to 
celestial latitude (3^ and longitude measured from the sun in the form of com- 

puter tabulations with superposed contours. The unweighted sample exhibits familiar 
features: concentration of most of the radiants to a broad band about the apex con- 
taining the "toroidal," "sun, " and "antisun" clusters, and a relatively small cluster 
at the apex itself. 

The atmospheric sample, even though smoothed, is comparatively noisy because 
of the high weights given a few slow meteors. Nonetheless, one sees that the sun 
and antisun clusters remain significant but that the remainder are more generally 
spread over the sky. 

In the space sample, the sun and antisun clusters have split into northern and 
southern divisions, of which only the northern divisions can be reliably observed 
at Havana. The total distribution now consists of two broad concentrations, which 
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RELATIVE NUMBER 


302-081 



Figure 5-1. 


Distribution of no-atmosphere velocities. The plotted numbers are 
fractions of the samples within a 2-km sec - * interval in velocity, or 
the fraction over 72 km sec - *. The observed, atmospheric, and space 
samples are explained in the text. 
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Figure 5-2. Distribution of observed radiants in latitude and longitude from the sun. 

The tabulated figures are the observed number within cells extending 5° 
in latitude and longitude, smoothed by averaging blocks of nine adjacent 
cells. The radiant density per square degree in a cell is equal to 0. 04 
sec (3 multiplied by the tabulated number. Contours have been drawn at 
4. 0 radiants per square degree (solid line), 1.2 (dashed), 0.4 (dotted), 
and 0. 12 (crosses). Contours were not drawn in southern latitudes, 
where the sample is obviously incomplete. The average density of 
observed radiants north of the ecliptic is 0.79 per square degree. 
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Figure 5-3. Distribution of radiants from the atmospheric sample, as in Figure 5 2. 

The weights have not been normalized, so that the contours represent 
only relative radiant densities. On the same scale, the average radiant 
density north of the ecliptic is 1. 40 per square degree. A few very 
high-weighted radiants have made this a rough distribution in spite of 
the smoothing. The obvious rectangular blocks of nine adjacent figures 
larger than their neighbors each represent one high-weighted radiant. 
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peak in these clusters. Moreover, the clusters have shifted away from the apex, and 
the relative numbers in the apex and antapex hemispheres are more nearly equal than 
any previously reported. 
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Figure 5-4. Distribution of radiants from the space sample, as in Figures 5-2 and 
5-3. On the scale of weights used, the average radiant density north of 
the ecliptic is 2. 86 per square degree. The smoothing process has arti- 
ficially lowered the bottom line of tabulated densities by one-third, but 
a real decrease also occurs at low latitudes. 


5. 4 Orbital Elements 


Figures 5-5 to 5-10 show the distributions of the elements a, 1/a, e, q, q', and 
i. The observed sample differs from previous observations primarily in the added 
importance of low-velocity meteors, which perforce have orbits with low to moderate 
eccentricity, low to moderate inclinations, and perihelia and aphelia not far from the 
earth's orbit. 
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RELATIVE NUMBER 



Figure 5-5. Distribution of semimajor axes a. The scale on the abscissa has been 
broken to permit showing details near 1 a.u. The plotted numbers are 
fractions of the samples per 1-a.u. interval of semimajor axis, or the 
fraction greater than 10. a.u. 



302 - 08 ) 



Figure 5-6. Distribution of inverse semimajor axes 1/a. The plotted numbers are 
the fractions of the samples within intervals of 0. 1 in 1/a. 

The atmospheric sample is dominated by near-circular or moderately eccentric 
orbits of low inclination. Other types are numerically insignificant. 

The space-sample inclinations are also low, but the peak of the distribution is near 
7 0 rather than 0° . Eccentricities, however, peak near 0. 5. These are in harmony, 
as they must be, with the radiant distribution and with a predominance of low geocentric 
velocities directed approximately radially to and from the sun. The general trend of 
the distribution of perihelion distance is a necessary consequence of observation at the 
earth and of a broad spread of eccentricities. However, a rise above the trend near 
0. 2 and a dip below 0. 1 may be indicative, respectively, of the perihelion of some 
significant source of the particles, or of particle removal by solar heating, or of both. 
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The general trend of the distribution of aphelion is also a consequence of observation at 
the earth, but a modest rise above the trend in the asteroid belt is very probably sig- 
nificant for the origin of many of these meteors (Section 7 has other evidence on the 
origin). 
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ECCENTRICITY e PERIHELION DISTANCE q 

Figures 5-7 and 5-8. Distributions of eccentricity e and perihelion distance q. The 

plotted numbers are the fractions of the samples within inter- 
vals of 0. 1 in e or q. 
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APHELION DISTANCE q' 


Figure 5-9. Distribution of aphelion distance q\ The scale has been broken to show 
more detail inside 4 a. u. ; this causes an apparent sharp change in the 
slopes that does not actually exist. The plotted numbers are the fractions 
of the samples per 1-a. u. interval of aphelion distance. 
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Figure 5-10. Distribution of inclinations i. The scale on the abscissa has been broken 
twice to permit a convenient representation. The plotted numbers are 
fractions of the sample per degree of inclination. 

Tables 5-1 to 5-3 show relative distributions of 1/a versus e, 1/a versus i, and 
e versus i, respectively, for the space sample. Each table was separately scaled to 
a maximum entry of 800. In Table 5-1, the upper blank area would contain orbits 
with perihelion outside the earth's orbit, and the lower blank area would contain orbits 
with aphelion inside the earth's orbit. Border cells of the triangular array are only 
about half observable from the earth, and scattered values outside the triangle are 
caused by the earth's orbital eccentricity. The concentration to the top edge of the 
array reflects the concentration of perihelion distance toward 1.0. Tables 5-2 and 
5-3 show small but definite positive correlations between l/a and i and between e and 
i, respectively. These await interpretation, as do many other prominent features 
and most details of the distributions. 
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Table 5-1. Relative distribution of l/a versus e for the space sample. 
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Table 5-2. Relative distribution of 1/a versus i for the space sample. 
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Table 5-3. Relative distribution of e versus i for the space sample 
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6. SPACE DENSITY 


6. 1 Observed Orbits 


The space density of meteors observed at Havana from 1962 to 1965, when the 
radar system was much less sensitive than during the synoptic year, was computed by 
Southworth (1967). The only other computation of space density (Briggs, 1962) treats 
an equilibrium distribution of photographic meteors under the Poynting-Robertson 
effect. The present report follows the methods of Southworth (1967), but it uses the 
far superior data from the synoptic year and reaches a very different result. 

The space density of the observed meteors is computed numerically in straightfor- 
ward fashion. Space in the solar system is divided into bins, with divisions at inter- 
vals of 0. 1 in log^r, of 0. 1 in sin (3, and of 18° in X, where r is the distance from the 
sun in astronomical units and where (3 and X are heliocentric latitude and longitude, 
respectively. The computer follows around each orbit and computes, to a very good 
approximation, the time spent in every bin traversed. Adding these times (with 
cosmic weights) and dividing by the bin volumes gives relative space densities. Abso- 
lute space densities then follow from the local space density derived from meteor influx 
data. 


6.2 Unobservable Orbits 


Only orbits with perihelion inside the earth's orbit and aphelion outside it can be 
observed from earth; clearly, we need to take the unobservable orbits into account. 

To that end, we extrapolate the element distribution observed at the earth, as follows: 
We assume that a and e are independent of the angular elements i, gj, and £2. Tables 
5-2 and 5-3 show that a and e are in fact only slightly dependent on i, while regression 
of the nodes under planetary perturbations must have randomized w and f2 . The two- 
dimensional relative distribution of l/a and e, shown in Table 5-1, was smoothed and 
then approximately fitted by 
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P(l/a, e) = BE , 




where 


E = 75- 131 (e-0.775) (|e-0.775|) °' 6 , > 


log 10 B = 


0. 2b 
3b 


if b £ 0 , 

if b > 0 , 


b= 1.35 - 1/a - e 


J 


( 6 - 1 ) 


Figure 6-1 shows the areas in the 1/a— e plane that contain orbits reaching 1 a. u. , 
reaching R a. u. , and reaching both. The correction factor F(R) to transform observed 
space densities at R to total space density was computed as the numerical integral 
of equation (6-1) in the region reaching R, divided by the integral in the common 
region. Table 6-1 is the result. 

Table 6-1. Correction factor for observed space densities. 
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3.24 


It is evident that the correction factor F must be unreliable when it is large, 
since large corrections are derived from the extrapolated existence of many unob- 
servable orbits. Although equation (6-1) was chosen as being a plausible form, as 
well as being an approximate fit to the observations, the bulk of the mass is in orbits 
with small e, for which the distribution of 1/a can be only poorly extrapolated. Thus, 
we should treat with reserve the corrected space densities beyond, say, 0. 4 and 
2 a.u., although the general trends beyond may still have significance. 
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I/O 

4 



0 I 

Figure 6-1. Orbits having 1/a and e in the horizontally shaded area can be observed 
at 1 a. u. from the sun; orbits in the vertically shaded area can be 
observed at R a.u. 

6. 3 Densities 


Table 6-2 shows decimal logarithms of space density in arbitrary units, in 
three two-dimensional arrays. The top section shows the distribution in r and sin p, 
averaged for all longitudes, and the bottom section shows the distribution in r and X, 
averaged for latitudes between -30° and +30°. Both these sections are still affected 
by the observational selection of meteors with nodes at the earth’s orbit and also by 
appreciable sampling noise, but we anticipate they may be of future interest. The 
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central section shows the distribution in r and sin p when the observed value of w for 
each orbit is replaced by a uniform distribution from 0 to 2tt. Regression of the nodes 
must in fact smear w nearly uniformly in the real interplanetary distribution; accordingly, 
we have simulated the process as our last correction for observational selection. 

Because the distribution is now symmetric about p = 0, this portion of the table includes 
only positive values of sin p. 


-22 -3 

Figure 6-2 shows the averaged density distribution, using 4 X 10 g cm for 

•f- T-> /-» nvNn rt/A rl rtvi n i n n r\ wi Cl 4- A J 4 v> On <-»4i *7 ’ I n niiwTrn n h n T.n V. n s~\ v\ 

cl j v_/ uvnu xi/y iiv/ui oiiv^ oux cxo ^ouiiaauvu xix i • xiiL> uui voo iiav^/ x^^»v>xx 

smoothed near r = 1 to correct bumps, which were doubtless caused by inexactness in 
the fit of equation (6-1) to the data. Of course, we could make cosmetic changes to 
(6-1), but we would make no real advance in reliability. It is not plausible that the 
earth should have any effect on space density that could be seen on the scale of 
Figure 6-2, especially at higher latitude. 


Figure 6-3 shows, on two scales, the contours of space density in a plane normal 
to the ecliptic. The less reliable contours are dashed. The meteors are concen- 
trated to a relatively thin layer about the ecliptic plane. There, the minimum density 
is 0. 7 a. u. from the sun, and the maximum, between 2 and 3 a.u. from the sun. 

6. 4 Comparisons with Other Data 

The distribution of space density just described was a considerable surprise to 
the authors. Nonetheless, it is in qualitative agreement with two recent results from 
totally different observing systems: (1) Roosen (1970) was unable to detect the earth's 
shadow on the gegenschein and deduced that the density of light- scattering particles 
must increase outward from the earth. (2) Pioneer 10 (Kinard and Soberman, 1972) 
detected an increase in particle density of larger particles on its way to the asteroid belt. 
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SOLAR DISTANCE (o.u.) 


Figure 6-2. Space density of meteors, corrected for unobservable orbits. Curves 
are plotted for 10 ranges of heliocentric latitude p and have been 
smoothed near 1 a.u. in dotted lines. 
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7. COLLISIONS IN SPACE 


7. 1 Significance of Collisions 

The drastic revision of meteor masses and orbital distributions communicated in 
earlier sections of this report led to a reappraisal of the importance of collisions 
between meteors in space compared to other loss processes, primarily the Poynting- 
Robertson effect. Formerly, it had appeared that collisions were important only for 
retrograde orbits, for large meteoroids, and for meteoroids near the sun. Now, after 
revision of the masses and therefore of the total flux by means of the new ionizing 
efficiency, it appears that collision is at least competitive with the Poynting-Robertson 
effect for ordinary orbits near the earth. 

The last previous estimate of radio meteor space densities (Southworth, 1967) had 
shown density decreasing monotonically outward from the sun, entirely in harmony 
with dominance of the Poynting-Robertson effect over all other loss mechanisms for 
meteors (except near the sun). The run of space density found in Section 6, however, 
can be maintained as a stationary distribution only if collisions are significant. 

7.2 Theoretical Collision Model 


A simple theoretical model for the distribution of space density when both collisions 
and the Poynting-Robertson effect are significant can be constructed as follows: Con- 
sider particles all of one size and in nearly circular orbits at small inclinations to the 
ecliptic. At some distance from the sun, the particles are supplied to the system at a 
uniform rate and then spiral in under the Poynting-Robertson effect, while preserving 
their respective orbital inclinations. Some of the particles at each distance from the 
sun undergo mutual collisions; we shall suppose that they are eliminated from the sys- 
tem through solar radiation pressure that blows their fragments away. The Poynting- 
Robertson effect reduces the radius r of each orbit at the rate 


(dr) 

Vdt / 


PR 


-5 X 10 11 
s 6r 


-1 

cm sec 


(7-1) 
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where t is time and s and 6 are particle radius and density in cgs units; the particles 
are taken to be spheres. If there were no collisions, the Poynting-Robertson effect 
would set up a density distribution inversely proportional to r (Southworth, 1967) 
because the orbital planes near the sun are more densely concentrated to the ecliptic. 
Thus, 



(7-2) 


Two particles collide if their centers approach closer than 2s, so that each 

o 

particle sweeps out a volume 4 tt s v dt in the time dt, where v is the mean relative 
velocity. Consequently, collisions diminish the space density at the rate 



4tt s 2 





where m= (4/3)ir so is the particle mass. 


(7-3) 


To determine v, consider first the mean relative velocity resulting from the 
orbital inclinations. Since these inclinations are small, the relative velocities are 
primarily normal to the ecliptic, and the mean absolute velocity component is 

v fl = 7.3X 10^ ir'^ 2 cm sec’ 1 , (7-4) 

where r is in centimeters and i is in radians. The mean relative velocity between 
particles is of the order of v , but depends on the distribution of inclinations. For 
this model, it will be sufficiently realistic to assume that inclinations are distributed 
from 0 to I so that the space density at the ecliptic plane is uniform for each interval 
of inclination. It can then be readily shown that the frequency of particles with 
inclination i is proportional to i itself, that the mean inclination is 

i = (2/3) I , (7-5) 
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and that the mean relative velocity normal to the ecliptic is 


v. = 

1 


MI V 

45 i n 


(7-6) 


Orbital eccentricity introduces relative velocity components parallel to the 
ecliptic plane, both radial from the sun and parallel to the circular velocity. The 
mean parallel component is approximately one-half the radial component. Since 
eccentricity also gives rise to collisions between particles at different mean distances 
from the sun, we must assume that the space density does not vary substantially 
between perihelion and aphelion of most of the particles. If we assume a con- 
tinuous distribution of eccentricities, it has not proved possible to obtain an analytic 
value for the mean relative velocity caused by eccentricity, but for a single eccen- 
tricity e we find the mean component parallel to the ecliptic plane to be 


v 


e 



(7-7) 


For a distribution of eccentricities, the numerical coefficient in (7-7) would be larger, 
because the faster particles collide more frequently and are therefore more highly 
weighted in the mean. When the mean velocity for a distribution of inclinations is com- 
pared with the mean for a single inclination, it appears that the coefficient should be 
increased by roughly a factor of 2. We therefore use 


3e 

v = — v 
e i n 


(7-8) 


where e is the mean eccentricity. Then for the total mean relative velocity, we have 


v= 7.3 X 10 12 (4.6 i 2 + 9 e 2 ) 1//2 r~ 1//2 


(7-9) 


Since there is a constant rate of particle supply, there is an equilibrium density 
distribution. As each group of particles supplied within a small interval of time 
spirals in toward the sun, the density within that group traces out the general run of 
density with respect to distance from the sun. Combining collisional losses (7-3) 
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with orbital shrinkage (7-1), and adding the effect (7-2) of the clustering of the orbital 
planes, we thus find that the space density obeys 


®pr + ®c®pr ' 


Substituting (7-1), (7-2), (7-3), and (7-9) gives 


d£= 

dr r 


+ A 


2 1/2 
P * 


9 


(7-10) 


(7-11) 


where 


A = 93 (i 2 + 2 e 2 ) 1 ^ 2 


(7-12) 


We now obtain the model density distribution in the ecliptic as the solution of 
(7-11): 



(7-13) 


which has a density minimum P m j n at distance r ^^ from the sun given by 


^min 


nun 


3 

2 


and 


r 

min 



, 2 1/2 
+ 3 r 0 




(7-141 


(7-151 


where is the density at a given distance r^ from the sun. Figure 7-1 shows the 
theoretical distribution (7-13). 
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Figure 7-1. Radial distribution of space density. The model distribution given by 

equation (7-13) is drawn as the smooth curve labeled "Simplified Theory" 
referred to the top and left-hand scales. The asymptotes are dashed. 

The top curve of Figure 6-2, labeled "Meteors" referred to the bottom and 
right-hand scales, has been positioned to make an approximate fit to the 
model. 


The model distribution (7-13) has a vertical asymptote at r = 2.25 r mj i) , corres- 
ponding to dominance of collisions over the Poynting-Robertson effect; but this part 
of the solution must be neglected because it violates the assumption that space density 
does not vary steeply with distance from the sun. The model distribution is also 
asymptotic at small r to 
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3 pr = p . r . 
r r min min 


(7-16) 


which corresponds to dominance of the Poynting-Robertson effect over collisions. 

We note that, in this simple model, if collisions eliminate most of the particles, 
then the asymptotic density at small r is related to the distance of the particle supply 
from the sun. Substitution of (7-14) and (7-15) into (7-16) shows that if 


p r « p r 
■asy asy r sup sup 


(7-17) 


then 


r = (2A p r ) 
sup ' K asy asy' 


-2 


(7-18) 


where subscripts sup and asy, respectively, denote the point of supply and a point 

on the asymptote. If we relax the condition (7-17), then the right-hand side of (7-18) 

is an upper bound to r 

* sup 


7. 3 Comparison with Observation 


For an order-of-magnitude comparison of the model distribution with observation, 
we can take Whipple's (1967) influx distribution and a mean geocentric velocity of 
9 km sec”^ (derived from our observed mean v^ of 14. 5 km sec *) to estimate a total 
space density of small particles (smaller than meteorites) of 


p 0 » 4X 10" 22 g cm" 3 (7-19) 

at the earth's orbit; thus, 

r Q = 1. 5X 10 13 cm . (7-20) 

Our observed inclination and eccentricity distributions give 

A *66 , (7-21) 
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and substitution into (7-15) gives 


r . ~ 1. 2 X 10 13 cm 
mm 


* 0.8 a. u. 


(7-22) 


Figure 7-1 also shows the space density on the ecliptic as derived in Section 6. 3, 
positioned to match at p Q and r m - n . We concentrate our attention on the interval from 
0. 5 to 1. 5 a. u. (and we must ignore the derived density at 0. 9 a.u. , as explained in 
Section 6). Within these limits, the fit between model and observation is better than 
might have been expected. In particular, the match along the r axis, given by (7-22), 
shows that we have been lucky in our estimates of A and p^. Nonetheless, observa- 
tion seems to support the approximate validity of the model in this interval. 

7.4 Sources of Meteors 

Comparison of the observed run of density in Figure 7-1 leads directly (if we 
accept the assumptions of the model and the reliability of the data analysis) to con- 
clusions concerning where the particles are supplied. No significant number of 
particles enter the system in the range where the model distribution fits the observa- 
tions - that is, from 0. 5 to 1. 5 a. u. Most of the particles enter between 1. 5 and 
roughly 4 a.u. - that is, in the region of the asteroid belt. However, this is also 
the region most densely populated by short -period comets and also by particles that 
are captured by Jupiter from long-period orbits, so that we cannot identify the source 
of this supply in an easy fashion. If we take the "observed" curve in Figure 7-1 at 
face value, there is a further supply of particles in the neighborhood of 0. 2 a. u. , 
which tempts us to identify them with particles ejected from bright long-period comets 
near perihelion. However, because of the limitation on the observations, it would not 
be safe to place any reliance on such a source at this time from this set of data. For 
the same reason, we cannot safely say anything about particle sources or their absence 
beyond 4 a. u. 

7. 5 Further Research 

Our recognition that collisions are very important in the evolution of meteor 
orbits requires us to recast the theory of orbital distribution. This is both an important 
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opportunity and a considerable task. The theory will need to include particle density, 
breaking strength, and fragment size distribution, as well as space density and flux 
distribution, because these are the parameters that determine the loss of most 
particles from the system and also the gain of collisional fragments to the system. 
Considerable computation will doubtless be necessary, but careful formulation is even 
more important. One result to be hoped is that the physical properties of individual 
particles can to some extent be determined from the other parameters, such as orbital 
distribution, that we now know to be linked with them. 

The simplified model in Section 7.2 neglects differences in collision rate that 
depend on particle size or orbital inclination or eccentricity. One effect of these 
differences is a difference in the average age of particles in the system, depending 
on size and orbit; large particles and those in highly inclined or eccentric orbits are 
much younger. Now that a meaningful space density distribution is available, it may 
well be possible to find collisional ages for meteor streams. 

At least three kinds of observational data are at hand for comparison with a 
theory including collisions, so that we have reason to expect that the comparisons 
can in fact tell us some things we should like to know, such as particle properties. 

First, the space density distribution is approximately as concentrated to the ecliptic 
plane as are the short-period comets, even though the long-period comets are the 
known source of both a substantial proportion of the observed large (photographic) 
and small (radar) meteors and almost all the observed dust input. Both Jupiter per- 
turbations (which are independent of the physical properties of the particles) and 
collisions and the Poynting-Robertson effect (which depend on the physical properties) 
are important to the present distribution. Second, there is vast observational 
literature on the mass distributions of the whole cloud and of individual meteor streams. 
The mass distributions of new streams are important input data to a general study of 
collisions; the overall mass distribution is a function of the physical properties of 
the particles. Third, the observed distribution of masses and orbits within individual 
meteor streams (including the 200 new streams announced in this report) makes it 
possible to compare dynamical and collisional ages of the streams. 
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Theoretical study of collision processes is, in fact, an essential step in improving 
our knowledge of interplanetary density distributions, because the emprical extrapola- 
tion of the observed orbital distribution in Section 6 needs theory if it is to be im- 
proved. At the same time, we would extend our knowledge of both general space 
distributions and the physical properties of individual meteors. (In particular, we might 
find that mass distributions depend on mean distance from the sun and perhaps from 
the ecliptic plane. ) We should also expect new results on the history of the sources of 
the interplanetary cloud. 
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8. COMPUTERIZED SEARCH FOR RADIO METEOR STREAMS 


8. I Two-Phase Computerized Stream Search in the Synoptic-Year Sample 

The stream search performed by us on the sample of 19, 303 radio meteor orbits 
from the years 1961 to 1965 was based on the statistical model of meteor streams 
(Sekanina, 1970a; referred to hereafter as Paper I). The model provides parameters 
of the distributions of meteor orbits in a stream and in the surrounding sporadic back- 
ground in terms of the Southworth- Hawkins (1963) D-criterion, which tests the 
similarity of meteor orbits in terms of the differences in their Keplerian elements. 

The search procedure starts from an initial orbit, which is assumed to be known, 
and iterates the elements of the stream's mean orbit by weighting them with the use of 
the elements of individual meteors until the orbit converges. The method is very 
convenient for searching for members (in a sample) of known meteor streams (see 
Sekanina, 1970b; Paper II) or for meteor streams associated with objects of known 
orbits (see Sekanina, 1973; Paper IH). Searching for new streams, for which of 
course no orbital information is available in advance, becomes more difficult as the 
size of a meteor sample to search grows. 

Theoretically, one should have to take the orbit of every single meteor in the 
sample as the initial set of elements, which for extensive samples is practically 
impossible, even with the aid of a powerful computer, because of the enormous 
number of calculations associated with the particular search procedure. 

From the synoptic-year data, 19, 698 radio meteor orbits have been successfully 
processed. This is far more than the number for which the orbit-to- orbit approach 
would be reasonable. Fortunately, it proved possible to use the original Southworth- 
Hawkins (1963) search program in the place of a preparatory search to supply the 
required initial orbits for further analysis by the statistical -model search program. 
Prnctical calculations have confirmed that the two computer programs complement 
each other excellently and have shown that the two-phase technique is the most effec- 
tive stream-search method in existence, particularly for a huge sample of orbital data. 
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We also used the orbits of previously known meteor streams and other objects as 
additional initial orbits, in part to check the completeness of the list of initial stream 
orbits supplied by the Southworth-Hawkins search program. 

8. 2 Phase I; Major Source of Initial Orbits 

With the synoptic-year sample of 19,698 radio meteors, the Southworth-Hawkins 
search program would require a computer memory capacity several times larger than 
that of a CDC 6400, on which this phase was carried out. It was therefore necessary 
to divide the whole sample into smaller sets and to search for streams in each of 
them. The sets had to overlap each other significantly to ensure that no stream was 
lost because of the dividing lines. This brought the estimated number of necessary 
subsets to at least 10. The longitude of perihelion and the inclination were taken as 
the dividing parameters, although the perihelion distance and the eccentricity were 
another possible pair. Only the nodal longitude and the argument of perihelion were 
felt to be less convenient, since for certain types of orbits, fairly large spreads in 
the two elements have only slight effects on the magnitude of the D-criterion. 

Table 8-1 shows the division of the sample into 12 sections by inclination and 

longitude of perihelion. For each section, we list the number of streams detected 

by the Southworth-Hawkins program at two different rejection levels, D . The 

s 

rejection levels were established in accordance with Lindblad's (1971) paper. The 

average number of meteors per section was about 3800, which, with Lindblad's 

formula (2), gives D = 0. 10. We also made an independent search with D = 0. 08 
s s 

in an effort to increase our ability to separate streams. For low-inclination 
sections, the rejection levels were decreased to 0. 08 and then 0. 06, because the 
major, broad streams are well known (we used the known orbits of these streams as 
the initial orbits to collect their members) and because we more urgently wished to 
separate possible minor showers. 

Table 8-1 indicates that 90 to 98% of the detected streams had no more than 
five members. We considered these streams as fortuitous groupings of nonrelated 
meteors and discarded them. By so doing, we reduced the total numbers of streams 
from 6019 to 502 for the higher rejection levels and from 4763 to 167 for the lower. 
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The remaining 669 streams went through a second screening: rejection of streams 
that were multiply listed because the 12 sections overlapped (Table 8-1) and because 
two rejection levels were used. We discarded all the multiple entries that differed 
from each other by less than 0. 06 (but keeping the limit somewhat flexible in individual 
cases). We also discarded entries that indicated a clear identity with a stream 
detected in Paper II or HI. We refrained from being too harsh in reducing the number 
of stream entries because multiple entries still surviving in the list would be rejected 
anyway in the following phases of screening. The present phase brought the number 
of streams down to 285. 

8. 3 Phase 1: Other Sources of Initial Orbits 


Besides the initial orbits provided by the Southworth-Hawkins program, we used a 
number of orbits of previously known streams and of potentially related objects as 
additional initial orbits for Phase II. Specifically, the entries are as follows: 

A. Orbits of the radio meteor streams detected in the 1961—1965 set of data, 
as published in Papers n and HI. 

B. Predicted orbits of twin showers (i. e. , daytime branches) of low-inclination 
streams of Papers II and III (see below for details). 

C. Orbits of periodic comets with perihelion distances q < 1. 2 a.u. , approaching 
the earth's orbit within a few tenths of 1 a. u. 

D. Orbits of earth- and Mars-crossing asteroids (with q & 1. 2 a.u. ). 

E. Orbits of meteorites and bright fireballs. 

Since verification of the existence of meteor streams detected in the 1961-1965 
list of orbits is one of the primary aims of this investigation, orbits A permit this to 
be done most easily. Also, the problem of stream associations with periodic comets, 
asteroids, and related bodies can be solved most directly by applying orbits C, D, 
and E, respectively. 

Orbits B deserve special attention. As pointed out by Whipple (1940) in the case 
of the Taurids, a broad meteor stream moving in a low- inclination orbit meets the 
earth at two points: on its way to as well as on its way from the sun. It is observed 
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as two apparently independent but actually closely related showers. In Paper III, we 
suggested that such a pair be called the twin showers. If the motion of the stream is 
direct, one of the twin showers (preperihelion branch) must be nighttime, and the other, 
daytime. In Table XIV of Paper III, we listed seven pairs of twin showers, of which 
three had been established a long time ago and the remaining four were detected in 
the 1961—1965 collection of radio orbits. 


In an effort to increase the chance of discovering more twin showers to the 
previously known nighttime low-inclination showers, we predicted approximate orbits 
of the former from the 1961-1965 orbits of the latter, using a criterion based on 
Figure 8-1. 

A stream moving in the plane of the ecliptic intersects the earth' s orbit on the 
way to the sun at point £2, and on the way from the sun at point £2*. 

For some members of the stream, £2 is actually the ascending node and w is the 
argument of perihelion; these meteors could cause a nighttime shower. For some 
other members of the stream, £2* is the ascending node and co* is the argument of 
perihelion; they could produce a daytime shower. Obviously, the two elements of 
the twin showers are related to each other by 


00^= 360° — co , 


£2 + — £2 + 2gj , 


( 8 - 1 ) 


while for the other three elements, we have i* ® i « 0, q* ~ q, e* » e. The predicted 
elements of the twin showers served as additional initial orbits in the Phase n 
calculations. 


8.4 Phase II: The Weighting System 

The technique developed in Paper I and previously applied in Papers n and III 
was used in Phase n with only one minor change. 


< 
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Figure 8-1. The orbits of twin showers (schematic). Here, to and S2 are the argument 
of perihelion and the longitude of ascending node, respectively, of the 
nighttime (preperihelion) shower, and to* and £2*, the same elements of 
the daytime (postperihelion) shower. 
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The search in the synoptic year sample for meteors belonging to 40 streams detected 
in the 1961—1965 collection was carried out with four different weighting systems. The 
system used for the 1961 — 1965 set, 


w = 1 - D/0. 2 for D < 0. 2 , 

w = 0 for D > 0. 2 , 

was generalized to 

w= (1 - D/D 0 ) h for D < D 0 , 

w = 0 for D > Dq , 


( 8 - 2 ) 


(8-3) 


and, with use of the same initial orbits, the mean orbits of the streams were cal- 
culated for the following combinations of (D^, h): 0. 20, 1; 0. 20, 2; 0. 25, 1; 0. 25,2. 

It was found that the four resulting sets of elements were in very good agreement 
for a number of streams and in reasonable agreement for most streams. 

There were, however, several instances where the mean orbits appeared to 
depend more substantially on the weighting system, and these were the key cases for 
deciding which of the four (D Q , h) pairs to prefer in searching for new streams. 

To illustrate the following arguments, we list a few examples in Table 8-2. For 
high-inclination streams, the mean elements from different weighting systems com- 
pare with each other quite favorably (December Draconids and X Draconids). At 
moderate and low inclinations, a high degree of consistency exists only for concen- 
trated streams (Geminids, Southern 6 Aquarids, and Monocerids). 

The system with = 0. 2 and h = 2 tends occasionally to give somewhat erratic 
results for very extended streams (Piscids, Southern Arietids, and to some extent £ 
Perseids), apparently because it gives too much weight to orbits with very low D's and 
too little weight to those with relatively large D's. 

The system with D^ = 0.25 and h = 1, on the other hand, seems to give too much 
weight to orbits with large D's, which may result in suppressing dispersed, 
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Table 8-2. Comparison of mean orbits of several streams calculated with different 
weighting systems (Dq,1i). 


Equinox 1950. 0 


Stream 


Southern 
6 Aquarids 


\ Draconids 


Triangulids 


Monocerids 


December 

Draconids 


D o 

h 

CJ 

Q 

i 

(a.u.) 

0. 20 

1 

60? 1 

81?2 

6? 6 

0.361 

0.20 

2 

62.9 

76.8 

4. 8 

0.385 

0.25 

1 

60. 1 

82.7 

7.6 

0. 364 

0.25 

2 

60.5 

80. 8 

6. 5 

0.365 

0.20 

1 

232. 1 

282.9 

0.7 

0.269 

0. 20 

2 

53.0 

102.5 

0.6 

0.280 

0.25 

1 

59.9 

83.0 

7.6 

0.363 

0.25 

2 

52.3 

102.0 

0.3 

0.274 

0.20 

1 

155.3 

305.7 

28.3 

0.070 

0.20 

2 

155.4 

305.4 

28.6 

0.067 

0.25 

1 

155.3 

306.2 

27.7 

0.071 

0.25 

2 

155.4 

305.7 

28.2 

0.069 

0.20 

1 

306.4 

172.5 

3. 8 

0.311 

0.20 

2 

297.0 

166.3 

5.6 

0.368 

0.25 

1 

306.7 

174.8 

2.0 

0.306 

0.25 

2 

306.3 

172.8 

3.5 

0.311 

0.20 

1 

155.9 

195.6 

66.4 

0. 967 

0.20 

2 

157.6 

195.9 

66.0 

0.973 

0.25 

1 

154.5 

195.5 

67.3 

0.962 

0.25 

2 

156.0 

195.7 

66.6 

0.967 

0.20 

1 

122. 1 

18.7 

2. 9 

0. 337 

0.20 

2 

133.2 

359.3 

2. 8 

0.239 

0.25 

1 

124.6 

14.2 

1.9 

0.315 

0.25 

2 

122.5 

17.8 

2.9 

0.333 

0.20 

1 

270.2 

218.4 

18.6 

0.603 

0.20 

2 

275.0 

215.6 

16.8 

0. 563 

0.25 

1 

264.2 

219.6 

18.7 

0.649 

0.25 

2 

269.6 

218.6 

18.3 

0.606 

0.20 

1 

135.4 

71.9 

22.3 

0. 155 

0.20 

2 

136.6 

71.4 

22. 1 

0. 147 

0.25 

1 

135.2 

72. r 

22.6 

0. 160 

0.25 

2 

135.8 

71. 8 

22.3 

0. 153 

0.20 

1 

185. 1 

253.7 

42. 1 

0. 981 

0.20 

2 

184.5 

253.6 

42. 8 

0.983 

0.25 

1 

183.3 

255.4 

41.8 

0. 981 

0.25 

2 

184.6 

254.2 

42.2 

0. 981 

0.20 

1 

325.3 

261.3 

23.4 

0. 139 

0.20 

2 

325.3 

261. 9 

23. 1 

0. 139 

0.25 

1 

325.2 

260.6 

23.4 

0. 140 

0.25 

2 

325.2 

261.4 

23.2 

0. 139 


n (1950) 
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low-population streams and lead relatively often to spurious convergences. Sur- 
prisingly, our calculations showed that this weighting system suppressed the (3 Taurid 
stream in favor of the £, Perseids. 

The two remaining weighting systems used give practically identical results in 
almost all cases. One of the very few exceptions was again the (3 Taurid stream, 
where the (0.2, 1) system gave the opposite node to that given by the other systems. 
The (0. 25, 2) system appeared to be slightly preferable for more dispersed streams 
because of being somewhat faster in their delineating, and just as good as the (0. 2, 1) 
one for more concentrated streams. Also, the quadratic form of the former approxi- 
mates a Gaussian-type curve, which is the postulated stream -membership probability 
function, much better than does the plain linear form of the (0. 2, 1) system. As a 
result, we found the (0. 25, 2) weighting system most convenient and applied it 
throughout. 



9. RESULTS OF THE STREAM SEARCH 
9. 1 List of the Detected Radio Streams 

Detection of any particular stream in different samples depends on the degree 
of compatibility of the corresponding observing schedules. The observing schedule 
of the synoptic year is given in Table 9-1, for comparison with the 1961- 1965 observ- 
ing schedule (see Table I of Paper I) . 

The orbits of possible streams provided by the Southworth-Hawkins search 
program, plus the orbits of potential parent or related bodies, of the streams detected 
in the 1961-1965 sample, and of suspected twin showers, produced a total of 473 
entries for Phase II of our search and resulted in detecting in the synoptic-year 
sample 256 streams, 200 of which were not found in the previous radio meteor sample. 

The orbital elements and related parameters of the streams detected in the 
synoptic-year sample are listed in Table 9-2. The individual columns give the follow- 
ing information (all angular data are referred to the standard equinox 1950. 0): 

Column 1. Designation of the stream. 

Column 2. Argument of perihelion gj and its mean error (degrees). 

Column 3. Longitude of ascending node £2 and its mean error (degrees). 

Column 4. Inclination i and its mean error (degrees). 

Column 5. Perihelion distance q and its mean error (a. u. ). 

Column 6. Eccentricity e and its mean error. 

Column 7. Semimajor axis a (a.u. ) and revolution period P (years). 

Column 8. Longitude of perihelion tt = oj + £2 (degrees) and the date of passage 
through the node (UT and days in 1950). 

Column 9. No-atmosphere velocity and its mean error (km sec *). 
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Column 10. Right ascension a of the corrected mean radiant from individual 

R 

meteor radiants and its mean error (degrees). 

Column 11. Declination 6^ of the corrected mean radiant from individual meteor 
radiants and its mean error (degrees). 

Column 12. Right ascension and sun-oriented celestial longitude X^-Xq of 
the corrected mean radiant from the stream's mean orbit (degrees). 

Column 13. Declination 6^ and celestial latitude (3^ of the corrected mean radiant 
from the stream's mean orbit (degrees). 

Column 14. Geocentric V^, and heliocentric Vjj velocity at the node (km sec - '*'). 

Column 15. The node (A — ascending; D — descending), and the type of stream 
(C - circumpolar, with the mean radiant permanently above the horizon at Havana; 

D — daytime; M — mixed, with the mean radiant never more than 10° below the 
horizon; N - nighttime). 

Column 16. Height at maximum ionization h and its mean error (km). 

& max v ' 

Column 17. Number of meteors in the sample with D < 0. 25 with respect to the 
mean orbit, with an asterisk to indicate that the stream was detected in the 1961 — 1965 
sample. 

9. 2 Identification of Streams Detected in the 1961— 1965 Sample 

The D-test was used to identify in the synoptic-year sample the streams previously 
found in the 1961— 1965 sample. With an identification threshold of D. = 0. 2, we were 
able to identify 55 of the 83 streams of the 1961—1965 sample; 1 more stream was 
probably detected, though with D. = 0. 26; and 27 streams were not found. The doubt- 
ful case is the Quadrantid shower. With the radar network idle during the period of 
the shower's maximum activity, the search program detected a minor stream with 
the mean orbit similar to that of the Quadrantids, but with the node shifted by 13° for- 
ward, thus matching the nearest operational period of the radar. Of the 27 streams 
apparently absent in the synoptic-year sample, positively three (January Scutids, 
Andromedids, and December Ursids) and probably two more (p Geminids and November 
Cepheids) could not be detected, because of gaps in the operational schedule in the 
critical periods. In 5 cases of the 22 remaining unidentified streams (Equuleids, 
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Table 9-1. Observing schedule of the Havana meteor radar during the synoptic year. 


1968 

December 2—6 

1 A 
-l*x 

December 16—20 

1969 

January 13-17 
January 27 — 31 
February 10- 14 
February 24 
February 26-28 
March 10—14 
March 17—19 
March 23-28 
April 7 — 9 
April 11-12 
April 17 
April 20—21 
April 24—25 
May 5—9 
May 19—23 
June 2 — 6 
June 11—12 


498 

2.53 

Old 

1 AO 

f-t -L^j 

X# vu 

269 

1.37 


808 

4. 10 

618 

3. 14 

496 

2.52 

75 

0.38 

89 

0.45 

557 

2.83 

25 

0. 13 

404 

2.05 

497 

2.52 

39 

0.20 

109 

0.55 

71 

0.36 

133 

0.67 

490 

2.49 

573 

2.91 

694 

3.52 

3 

0.01 


1969 (cont. ) 

June 16-21 

Tii«a O C\ TiiIt* a 

U Ull^ WU U Uij -X 

July 14—19 
July 21 
July 25—26 
July 28 -August 1 
August 11—15 
August 17— 18 
August 25—30 
September 7 
September 9—12 
September 18 
September 22—26 
October 6—10 
October 14—15 
October 20—24 
November 3—8 
November 14 
November 16 
December 12—14 

Total 


701 

3.56 

rv7 n 

sJ 1 1 

A QG 

"X* i/V 

1059 

5.38 

3 

0.01 

11 

0. 06 

1272 

6.46 

1373 

6. 97 

175 

0.89 

892 

4. 53 

219 

1. 11 

1239 

6.29 

49 

0.25 

1210 

6. 14 

1369 

6.95 

44 

0.22 

1117 

5. 67 

1124 

5.70 

8 

0.04 

39 

0.20 

157 

0. 80 

19, 698 

100. 00 


Date (CST) 


Number of 
successfully 
reduced meteors 


Date (CST) 


Number of 
successfully 
reduced meteors (%) 


59 




Table 9-2. List of 256 streams detected among 19,698 radio meteors of the synoptic-year sample. 
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X Piscids, April Draconids, k Cygnids, and (3 Cetids), the search revealed streams 
that must have been classified as new, but that differed significantly from the corres- 
ponding 1961—1965 streams in only one element other than the nodal longitude. 

Of the remaining 17 streams of the 1961—1965 sample missing in the synoptic- 
year sample, 11 were found to be "absorbed" by one of the 55 identified streams 
(X Capricornids, a Leonids, May Librids, May Herculids, v Librids, Taurids- 
Perseids, y Aquarids, a Cygnids, 6 Cygnids, X Draconids, and October Camelopar- 
dalids); 2 by the sporadic background (June Draconids and Andromedids); 1 by a 
new stream (ir Cepheids); and in 3 cases, no or only a few possibly associated meteors 
were found (Lyrids, rj Aquarids, and Halleyids). 

9. 3 New Streams; Toroidal Group, Short- and Long-Period Streams, Twin Showers , 

and the "Cyclid" System 

A great number of the new streams belong to the so-called toroidal group (Hawkins, 
1963). They have direct but high-inclination orbits with low to moderate eccentricities. 
For example, we find 10 streams with i > 70°, of which all have q > 0. 6 a.u. , but only 
3 have semimajor axes significantly larger than 1 a. u. Their radiants lie at high 
northern latitudes, mostly in the constellations of Draco, Ursa Major, Camelopardalis, 
Cassiopeia, Cepheus, and Ursa Minor. The traditional system of stream designation 
used in Paper III proved insufficient to provide names for all the new streams in that 
area of the sky so heavily populated by radiants. We therefore expanded the previous 
system of designation by adding capital letters before the name of the radiant's con- 
stellation wherever the original system proved inadequate. Also, streams coming 
from the constellation of Ursa Minor have been called Umids to distinguish them from 
streams from Ursa Major, which are called Ursids. 

There are only two retrograde streams in the sample, the Orionids and the 
Perseids. The latter were not classified as a detected stream in the 1961— 1965 
sample, since only two probable members were found there. 

A total of 26 streams have their perihelia less than 0. 2 a.u. from the sun, and 
9 streams less than 0. 1 a.u. from the sun, the record being 0. 045 a.u. for the 
November Orionids. 
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Nine streams have semimajor axes shorter than 0. 8 a.u. , and three shorter than 
0.7 a.u. Two of these, the J Camel opardalids and p Andromedids, actually cross the 
earth's orbit almost at their aphelion points (aphelion distances 1. 02 and 1. 06 a.u. , 
respectively). 

On the other hand, there are 12 streams in the sample with semimajor axes of 
3 a.u. or longer. Surprisingly, neither the Perseid stream, nor the Orionid, nor the 
Monocerid holds the record. Rather, a new, very rich stream with q = 0. 86 a. u. , 
inclination almost 50°, and revolution period of more than 30 years, appears to be 
the most extended, not only in the synoptic-year sample, but probably among all radio 
streams known so far. It is called the September Ursids, and it is associated with 
neither a known comet nor a photographic stream. 

We now know 14 pairs of twin showers with confidence and 4 more pairs possibly 
related. Two of the seven pairs listed in Paper III as detected in the 1961— 1965 sample 
were not found in the synoptic-year sample: the Orionid-Halleyid pair because of the 
absence of the Halleyid stream; and the cr Capricornid-X Capricornid pair because of 
the absence of the X Capricornid stream. The e Arietids replaced the X Piscids as the 
twin shower loosely associated with the Triangulids. The new list of detected twin 
showers is in Table 9-3; the D values for the uncertain cases are in parentheses. 

Nine of the new streams have inclinations less than 10° and eccentricities less 
than 0. 4. We find that there is a fairly close relationship among most of them in 
terms of D despite an apparent variety in their orbits (Table 9-4). A similar relation- 
ship detected by Southworth and Hawkins (1963) among orbits of this type of photo- 
graphic meteors led the two authors to call such a "stream" the Cyclids and interpret 
it as a result of a high probability that meteors in these orbits meet the earth. 

9.4 Distribution of Radiants 


The computer plots of radiants of the 19, 698 radio meteors of the synoptic -year 
sample are exhibited in Figures 9-1 and 9-2. Figure 9-1 is a plot of the ecliptical 
coordinates, longitude X versus latitude (3, in Hammer's equal-area projection. 
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Figure 9-2 is a plot, in the same projection, of the sun-oriented ecliptical longitude 
X - \q versus latitude p. Except for a few details, the plots are almost identical with the 
corresponding plots of the individual meteor radiants from the 1961—1965 sample (see 
Figures 4a and 4b of Paper III). 


Table 9-3. Twin showers detected during the synoptic year. 


Twin showers 

D 

i 

Preperihelion 

Postperihelion 

Southern Arietids 

£ Perseids 

0. 108 

5° 

Scorpiids-Sagittariids 

C apricornids-Sagittariids 

0. 119 

4° 

X Orionids 

k Aurigids 

0. 135 

3° 

Southern r) Virginids 

if> Virginids 

0. 136 

2° 

6 Leonids 

August Lyncids 

0. 144 

10° 

January Cancrids 

August Cancrids 

0. 148 

3° 

p Leonids 

t Leonids 

0. 150 

1° 

OphiucMds 

January Sagittariids 

0.150 

2° 

Piscids 

May Arietids 

0. 154 

3° 

Taurids 

P Taurids 

0. 157 

0° 

p Geminids 

£ Leonids 

0. 164 

2° 

y Piscids 

X Capricornids 

0. 182 

5° 

Librids 

£ Sagittariids 

0. 198 

3° 

a Taurids 

p Cancrids 

0.202 

1° 

6 Cancrids 

Y Leonids 

(0. 285) 

4° 

a Capricomids 

c Aquarids 

(0. 288) 

7° 

t, Aurigids 

March Andromedids 

(0. 320) 

12° 

Triangulids 

e Arietids 

(0.388) 

11° 


The mean radiants of the detected streams (both new and known from the previous 
sample) are plotted in Figures 9-3 and 9-4. They show much the same pattern as the 
plots of individual radiants. At least, there is much better agreement between 
Figures 9-2 and 9-4 than between Figures 4b and 5b of Paper in. This is of course 
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due to the fact that the present search was much more complete than the former. 
Although we have not yet been able to determine the population and dispersion coeffi- 
cients of the streams found in the synoptic-year sample, we can estimate that the 
total number of meteors in the detected streams is of the order of 6000, or approxi- 
mately 30%, as compared to 1400 or 7. 5% in the 1961— 1965 sample. 

Table 9-4. The Cyclid system of streams in the synoptic year (the numbers are D 
values in units of 0. 001). 


•a ® 

•n V. 


Stream 

CG 

1 

a 

k Geminids 

£ 

■8 

& 

2 

\ Aurigids 

a Aurigids 

June Aurigidi 

CG 

xs 

•l-l 

CG 

U 

P 

vJ> 

« Ursids 

p Ursids 

| Cygnids 

0 

188 

370 

254 

442 

493 

272 

266 

260 

k Geminids 

188 

0 

266 

93 

238 

267 

229 

216 

165 

March Virginids 

370 

266 

0 

324 

404 

387 

269 

244 

252 

X Aurigids 

254 

93 

324 

0 

164 

236 

322 

306 

252 

a Aurigids 

442 

238 

404 

164 

0 

149 

424 

466 

379 

June Aurigids 

493 

267 

387 

236 

149 

0 

449 

427 

309 

t, Ursids 

272 

229 

269 

322 

424 

449 

0 

87 

182 

« Ursids 

266 

216 

244 

306 

466 

427 

87 

0 

127 

p Ursids 

260 

165 

252 

252 

379 

309 

182 

127 

0 


9. 5 Comparison with Cook’s Working List of Meteor Streams 

Of Cook's (1972) working list, 16 streams were closely matched by radio streams 
from the synoptic-year sample (D H 0. 1), and 18 more streams were probably identi- 
fied (D SL 0.2). However, we were not able to separate the northern and the southern 
branches of the Taurids and the X Orionids. Because the radar system did not operate 
early in January, we found only a remote branch of the Quadrantids. The remaining 
21 streams of Cook's list, many ctf them either retrograde in motion or nearly para- 
bolic or ones in southern hemisphere, were not found in the synoptic-year sample. 
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Figure 9-3. Ecliptical coordinates X and p of the mean radiants of 256 radio meteor streams in Hammer's equal- 
area projection. Open circles stand for the Cyclid-type streams; circled points, for twin showers. 
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10. COMET- METEOR AND ASTEROID- METEOR ASSOCIATIONS 


10. 1 Associations with Periodic Comets 


Table 10-1 lists the detected associations between the meteor streams of the 
synoptic-year sample and the potential parent objects within D = 0. 3. 

We find that the following periodic comets are likely to have associated streams 
in the synoptic-year sample: Encke, Pons-Winnecke, Giacobini-Zinner, Halley, 
Schwassmann-Wachmann 3, Swift-Tuttle, and Mellish. Possible associations have 
also been found for some other comets (Honda - Mrko s - Paj du^ akova, Gale, Helfenzrieder, 
Blanpain, and Lexell), but these are suspected to be fortuitous, because of the comets' 
orbit instabilities in some cases or because of relatively loose relationships in others. 

10.2 Adonis Meteor Streams 


Among the asteroid-meteor associations, the most obvious one is again the case 
of Adonis. The cr Capricornids, a prominent stream in the 1961—1965 sample, were 
found to match the asteroid's orbit less satisfactorily in the synoptic-year sample, 
partly because of smaller perihelion distance and partly because of higher inclination. 
However, it seems increasingly likely that the Scorpiids-Sagittariids, a very broad 
stream undoubtedly related to Hoffmeister' s (1948) Scorpiid-Sagittariid complex, is 
associated with Adonis. We note that both the June nighttime stream and its twin 
shower, the Capricornids-Sagittariids appearing in January and February, are in 
terms of D within 0. 2 from the asteroid's orbit, and that a minor stream, the 
X Sagittariids, closely related to the Scorpiids-Sagittariids (actually its component), 
is within 0. 09 of the asteroid's orbit. Hoffmeister did notice a resemblance between 
the orbits of Adonis and his Scorpiid-Sagittariid system, but left open the possibility 
of a real association between both. This was the correct approach because his results 
were based on inaccurate visual observations; also, he found a discouraging difference 
of almost 25° between the longitudes of perihelion of the asteroid and of the middle 
of the meteor complex. Our radio meteor samples suggest that the difference is 


85 


Table 10-1. Possible comet-meteor, asteroid-meteor, and meteorite/fireball-meteor associations. 
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actually much less than that. In the synoptic-year sample, the Capricomids- 
Sagittariids indicate a difference of 13°, and the Scorpiids-Sagittariids, only 8°. 

10. 3 Possible Associations with Other Minor Planets of the Apollo and Amor Types 
and with Meteorites and Bright Fireballs 

Besides Adonis, there is a definite possibility of the existence of associations of 
meteor streams with the minor planets Apollo, Hermes, Toro, and 1950 DA, and to 
a lesser degree also with Icarus and, surprisingly, Geographos and even Eros. 

The last two must, however, be taken with much caution: The would-be associated 
streams belong to the group of Cyclids. We know that these orbits have high 
probability of encounter with the earth, and detection of any well-established associa- 
tions is most difficult if not impossible. Besides, the perihelion distance of Eros is 
1. 13 a.u. , so that detectable meteors could only be fairly loosely related to the 
asteroid anyway. 

Of the meteorites and bright fireballs, the Pribram meteorite and the Prairie 
Network fireball 40617 may be associated with some of the detected streams (see 
Table 10-1). No streams were found to move in the orbit of the Lost City meteorite, 
but the radar system was not operating at and around the node of that object. No 
streams seem to be associated with a few other fireballs, among them Prairie Network 
40503 and Mt. Riffler (in the latter case, the radar was, however, again idle). 
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11. HEIGHT -VELOCITY DIAGRAM 


11. 1 Introduction 


Ceplecha (1967, 1968) detected three discrete levels of beginning height on height- 
velocity plots of photographic meteors and interpreted them in terms of the composition 
and fragmentation of the meteoroids, thus suggesting that the bulk density of the latter 
is important. Cook (1970) applied Ceplecha's criterion to classify 25 photographic 
streams. 


Height-velocity computer plots have been used in this study to search for the dis- 
crete height levels in the synoptic-year sample. Since for radio meteors the height 
at maximum ionization is much easier to define than the beginning height, we used 
essentially the former. A few check plots with the beginning height have confirmed 
that it makes little difference which of the two heights is used, except that the begin- 
ning height is subject to a larger dispersion. 

11. 2 Height -Velocity Plots for Individual Radio Meteors 

The plot of 16, 322 radio meteors of the synoptic year with zenith distance of 
radiant less than 60° (to eliminate a possible effect of nearly ’’tangential" meteors) 
has failed to show the discrete levels detected by Ceplecha among the photographic 
meteors (Figure 11-1). Nothing has been gained by dividing the whole sample into 
groups according to the meteor luminosity and perihelion distance (Figures 11-2 to 
11-7). Dispersion seems to be large enough to mask completely the double-peak dis- 
tribution in height. However, we have been able to detect some of the features of the 
height-velocity plots known for the photographic meteors. 

First, there is a clear separation between the major concentration (at low Voo) 
and the minor concentration (at high V M ) in the velocity axis at ~ 48 km sec . 

The effect is obvious from all figures except Figure 11-4. In Ceplecha's terms, the 
gap separates the minor C 2 group from the major (plus A) group. The separation 
was originally recognized by Jacchia (1958). 


89 



HEIGHT AT MAXIMUM IONIZATION (km) 



9 ' 







HEIGHT AT MAXIMUM IONIZATION (km) 


70 


Figure 11 


80 


90 


100 


110 


10 


20 


30 


50 


80 


NO- ATMOSPHERE VELOCITY (km sec -1 ) 


-3- Height-velocity plot of 7835 radio meteors of the synoptic year with 
zenith distances less than 60° and "fainter" at maximum ionization 
than 11. 5 mag. 
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Figure 11-5. Height-velocity plot of 3008 radio meteors of the synoptic year with 
zenith distances less than 60° and perihelion distances between 0. 25 
and 0. 50 a.u. 
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Figure 11-6. Height -velocity plot of 3256 radio meteors of the synoptic year with 
zenith distances less than 60° and perihelion distances between 0. 50 
and 0.75 a.u. 



NO-ATMOSPHERE VELOCITY (km sec ') 


Figure 11-7. Height-velocity plot of 7181 radio meteors of the synoptic year with 
zenith distances less than 60° and perihelion distances larger than 
0. 75 a.u. 




Second, the vertical dispersion in the major concentration (V*, < 48 km sec *) is 
definitely larger than that in the minor concentration (V„o > 48 km sec *). From 
Ceplecha’s results, we actually expect this, because the dispersion in the major con- 
centration is a sum of the dispersions of the two groups (A + Cj) several kilometers 
apart, while the dispersion in the minor concentration represents that of the C 2 group 
alone. 

Third, the minor concentration is more pronounced (relative to the major concen- 
tration) among the brighter meteors. Ceplecha found the same effect among the 
photographic meteors, although the magnitude scale was naturally shifted by about 
10 mag. 

Finally, Figure 11-4 (meteors with perihelion distances less than 0.25 a.u.) is 
the only one that shows but one cluster of meteors with the cutoff at ~ 60 km sec -1 
rather than 70 to 75 km sec - '*' (as the other figures do). Among the photographic meteors, 
the group with q < 0.25 a.u. also indicated only one concentration; Ceplecha classified 
these meteors as forming an intermediate height level, which he called the B group. 

11. 3 A Height- Velocity Plot for Radio Streams 

The last attempt to separate the discrete levels of radio meteor heights was a 
height-velocity plot for the meteor streams from the synoptic-year sample. Only 
streams with a mean error in height not exceeding ± 1 km were plotted. The idea 
was that with individual heights grouped into the stream heights, the dispersion within 
each of the two parallel levels A and C ^ may be lowered, so that the two levels show 
up. That did not, however, turn out to be the case, as can be seen from Figure 11-8. 

The major concentration remains structureless, despite a decrease in the overall dis- 
persion along the height axis, as compared to the individual height -density plots. 

An interesting point, possibly worth pursuing further in the future, is a clear indica- 
tion that the area of the main cluster with maximum heights 93 to 96 km at 25 km sec \ 
«95 to 98 km at 35 km sec -1 ) is not populated by ecliptical high-eccentricity streams, 
as one might expect by analogy to Ceplecha' s conclusions, but mostly by high- inclina- 
tion, low -eccentricity streams. A question to be answered is whether a possible 
selection effect exists, since these streams consist of meteors whose radiants are 
most favorably located for the main beam of the radar system. 
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Figure 11-8. Height-velocity plot of 145 radio meteor streams with reliably deter 
mined mean heights (mean error in height not exceeding ±1. 0 km). 

Another feature of Figure 11-8 is a drastic drop in the height dispersion at 
velocities below about 18 km sec \ At the other end of the velocity range, the C 0 
group has reduced to only two points: the Orionids and the Perseids. 
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12. SPACE DENSITY IN RADIO METEOR STREAMS 


12. 1 Mean Space Density in a Stream in Terms of the Sporadic Density 

In Papers EE and HI, we derived the two parameters of the statistical model of 
meteor streams — the population coefficient A and the dispersion coefficient a — for 
a number of streams from the 1961— 1965 sample, and we commented that the popula- 
tion coefficient gives essentially the flux ratio of the stream meteors relative to the 
sporadic background in a limited vector space (i. e. , the stream-to-background flux 
ratio of meteors from within a particular solid angle related to particular limits of 
the geocentric velocity). By varying the value of the D-test, we change the "volume" 
of the vector space and thereby can study the variations in the stream-to-background 
flux ratio along a given section of the earth's orbit. By extending the vector space 
to cover all directions and all velocities, we could derive the stream-to-background 
space -concentration ratio. 

The parameters of the statistical model are derived from the cumulative D-distri- 
bution of meteors, not of meteor masses. This makes no difference within any stream, 
because the total range of D is then fairly narrow. Selection effects caused by both 
the observing conditions and the physical processes are, however, of considerable 
importance in assessing the overall concentration of meteoric matter in space, as 
shown elsewhere in this report. 

The mean space density in a meteor stream, in units of the mean space density 
of the meteoric matter in the sporadic background, was derived by us with the use 
of the stream's population and dispersion coefficients and with the use of the overall 
distribution erf the orbital elements of meteor masses from the synoptic year, corrected 
for all selection effects (see Figures 5-5 to 5-10). In Paper I, we defined the outer 
limit Djj of a stream by a condition of equal population of the stream and background 
meteors between D = 0 and D = Djj. The expression for D^ and its graph are also 
given in Paper I. Papers n and III list D^ for 83 streams of the 1961—1965 sample. 
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For each stream, a Monte Carlo method has been applied to simulate the distribution 
of 10, 000 sporadic meteors of equal mass with orbital characteristics distributed as 
given in Figures 5-5 to 5-10. This distribution gives us the background noise for each 
stream in terms of the overall distribution of meteor masses. Specifically, for 
D = Djp we can immediately find the fraction of the sporadic meteor masses that falls 
within Djj from the mean orbit of the stream. From the definition of D„, this is also 
the fraction of the number of meteors (or meteor masses) of the stream that falls 
within Djp which covers practically all the stream members. Since the space density 
is proportional to the total of meteor masses swept up by the earth within a certain 
period of time, this method gives us the space density in the meteor stream within 
Djj in terms of the overall space density of the sporadic background. 

12. 2 Space-Density Gradient in a Stream 

In any meteor stream that can be distinguished from the sporadic background, 
there exists a measurable space-density gradient. The method described gives 
actually a lower limit for the space density in a stream, because Djj is much larger 
than the stream's characteristic dimension (~ oV2"). We do not know the profiles of 
most streams, because there were gaps in the observing schedule, but we can 
estimate the density gradient from the shape of the D-distribution curves in the 
following way. 

The number of meteors with D < D., Jj? s ( D.), swept by the earth during a pass 
through a stream is 


3g,(D.) = const • p(D.) W(D.) cosec c , (12-1) 

E 

where p(D.) is the mean space density in the stream for D 2 = D.; W(D.) is the charac- 
teristic width of the stream for the same range in D; e is the elongation of the 
corrected radiant from the apex of the earth's motion; and V„ and V„ are the comet's 
and earth's heliocentric velocities, respectively. For p independent of D, we would 
have from equation (12-1), ~ W. From Paper I, we know that the number of 

meteors would then vary as dP (where (3=3. 8), so that W ~ D* 3 . In the case of a 
variable space density in the stream, equation (12-1) can be written thus: 
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Jg(D.) = const • p(D.) D? . 


( 12 - 2 ) 


It is convenient to replace D. as the argument by E. = D./(oV5~) and to write, with the 
use of equation (6) of Paper I, the ratio of the mean space density in the stream 
within D. to that within c rtt : 


P< E i> -# S < E i> -3.8 

Pd) PM) i 

‘ '-U 



1. 13 E { e 



e7 3 * 8 


where 


erf (E) = 




m-3) 
\ / 


(12-4) 


The function p(E)/p(l) is plotted in Figure 12-1. It can be shown that for E « 1, this 
expression simplifies to 


p(E) = 
P(l) 


1.77 


E 


- 0.8 

> 


and for E » 1, to 


£1E) 

P(l) 


= 2. 34 E -3 * 8 . 


(12-5) 


(12-6) 


Also plotted in Figure 12-1 is the relative number of meteors in a stream as a function 
of E. We note that the statistical model predicts that only 0. 1% of the total mass of a 
stream move in orbits with D < 0. 15 a, so that the mean density within E = 0. 1 is 
perhaps a realistic estimate for the actual maximum central space density in meteor 
streams. 


12. 3 Numerical Results 


The results of the Monte Carlo simulation, described earlier in this section, 
and those of the subsequent space-density calculations are summarized in 
Table 12-1 for several of the more important streams of the 1961—1965 sample. 


101 


We conclude that mean space densities in meteor streams are small compared with 
the space density of the surrounding sporadic background, but that central densities in 
some major streams can significantly exceed the background density. 


301 - 052 



Figure 12-1. Space-density variation (solid curve, left-hand scale) in meteor streams 
versus E = D/(aV2") (<T is the dispersion coefficient) predicted from the 
statistical model, and a relative cumulative number of a stream's mem- 
bers (dashed curve, right-hand scale) as a function of E. 

-22 -3 

With an average overall space density of 4 X 10 g cm in the vicinity of the 

earth (see Section 7. 3), the results of Table 12-1 suggest that mean densities in 

-24 -3 

meteor streams are typically several times 10 g cm and that the central density 

-20 -3 

in some streams could be as high as almost 10 g cm . In meteoric storms, the 

-19 

density must be still at least 1 to 2 orders of magnitude higher, i. e. , some 10 to 

m-18 -3 

10 g cm . 
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Table 12-1. Mean and central space densities in radio meteor streams. 


Stream 

Space density in stream 

-22 -3 

(in units of sporadic density =* 4 X 10 g cm ) 

mean (at D ^ D^) 

central (at D < 0. 15 a) 

Sc orpiids -Sagittar iids 

0.06 

30 

Geminids 

0. 005 

i n 

JL V 

a Capricomids 

0.05 

8 

Piscids 

0.03 

3 

Monocerids 

0.01 

3 

Taurids 

0.02 

2 

Southern Arietlds 

0.02 

2 

July Draconids 

0.006 

2 

Triangulids 

0.02 

1 

o Cetids 

0.006 

1 

Southern 6 Aquarids 

0.003 

1 

Lyrids 

0. 0002 

0.4 

Arietids (Daytime) 

0.002 

0.3 

X Orionids 

0.003 

0.2 

£ Perseids 

0.002 

0.2 

(3 Taurids 

0.002 

0. 1 

a Leonids 

0.002 

0.1 

Southern Virginids 

0.001 

0.05 

Quadrantids 

0. 00002: 

0.01: 

Orionids 

0. 000002: 

0. 005: 


12.4 Space Density from Cometary Production Rates of Solid Material; Comparison 
of the Two Methods 

It is of interest to compare the above conclusions with independent evidence avail- 
able from the study of cometary type II tails. Finson and Probstein's (1968a) model 
of dust comets makes it possible to estimate the total mass of the solid material lost 
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by a comet during an approach to the sun. The results from the regular dust tails are 

unfortunately inconclusive for the particle sizes that produce radio meteors (« 0. 1 cm 

in diameter), since the method is based on photometric considerations and the effect 

from heavy particles is always masked by the much stronger light- scattering power of 

micron and submicron particles. Consequently, the size-distribution function, which 

is also determined by the Finson-Probstein model, is uncertain for large masses. 

However, the "spike" tail of comet Arend-Roland 1957 ID, studied by Finson and 

Probstein (1968b), appeared to be composed of only heavy particles (> 0. 01 cm in 

diameter). The two authors showed that the spike had been formed as a result of an 

outburst about 2 months before the comet reached the perihelion point and that the 

13 14 

total mass ejected during the outburst had been of the order of 10 to 10 g. It 
seems that this is the only quantitative estimate to date of the emission rate of dust 
particles approaching the radio meteor range of masses. 

With a characteristic dispersion in their velocity distribution along the orbit of 

some 10 m sec the particles ejected from a short -period comet in an isolated 

outburst would form, after one revolution around the sun, a filament some 0.2 a.u. long. 

With a width of * 0. 001 a.u. (corresponding to the duration of less than 0. 1 day, 

32 3 

typical for meteoric storms), the volume of the filament is * 5X10 cm , and with 
the mass of « 10 ^ g, its space density comes out to exceed 10 ^ g cm in total 
agreement with the estimate given in the previous section. 

To estimate the space density in old streams dispersed along the whole arc of 

orbit around the sun, we can adopt the orbital length » 15 a.u. and the width = 0.5 a.u., 

40 3 

so that the volume of the stream is ~ 10 4U cm . The total mass involved can be 

3 

inferred either from the order-of-magnitude estimate of the age (= 10 revolutions), 
assuming that the comet is ejecting meteoroids at a more or less constant rate, or 
from the total mass of the parent comet (if the latter is presumed defunct), assuming 
that a significant fraction of the comet disintegrated into the stream. Both estimates 

lfi 17 _9 o 

give fairly consistently a mass of « 10 to 10 g and a space density ~ 10 to 
-24 -3 

10 gem, again consistent with the results of Section 12. 3. 
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